Intrinsic giant Stark effect of boron-carbon-nitride nanoribbons with zigzag edges 
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Electronic properties of zigzag boron-carbon-nitride (BCN) nanoribbons, where the outermost 
C atoms on the edges of graphene nanoribbons are replaced by B or N atoms, are theoretically 
studied using the first-principles calculations. We show that BCN nanoribbons are metallic, since 
several bands cross the Fermi level. For BCN nanoribbons in a rich H2 environment, the so-called 
nearly free electron state appears just above the Fermi level because of the intrinsic giant Stark 
effect due to the internal electric field of a transverse dipole moment. The position of the nearly 
free electron state can be controlled by applying an electric field parallel to the dipole moment. The 
hydrogenation of the nitrogen atom is necessary for the appearance of the giant Stark effect in BCN 
nanoribbons. We also discuss the effect of stacking order on the intrinsic giant Stark effect in bilayer 
BCN nanoribbons. 

PACS numbers: 73.20.At, 73.21. Cd, 73.22.-f 



Graphene nanoribbon (GNR) is a nanometer- wide 
strip cut from a graphene (Gr) sheet on which C atoms 
are arranged in a honeycomb lattice. The electronic and 
magnetic properties of GNR have been studied inten- 
sively by several groups [1-4 , which predicted the pres- 
ence of a peculiar electronic state at zigzag edges, which 
is called the edge state. GNR can be fabricated by the 
use of electron-beam lithography [5], the unzipping of 
carbon nanotubes (CNTs) |6l[7], and the use of bottom- 
up processes [8 . The presence of edge states of GNR 
were confirmed experimentally by Tao et al [9]. 

One of the most interesting properties of Gr-based 
nanomaterials, such as CNT and GNR, is the controlla- 
bility of energy band structures using an applied electric 
field, i.e., the giant Stark effect p!QHT5] . By applying a 
transverse electric field, the band gap of semiconducting 
CNTs can be closed jlOHIS] and the linear dispersion of 
metallic CNTs around Fermi level can be bent [M*, 15^ . 

Inspired by the rich physics and functionalities of 
CNT and GNR, researchers have studied nanotubes 
and nanoribbons made of hexagonal boron-nitride (BN) 
sheets, which are an inorganic analogue of Gr [16 . The 
giant Stark effects have been predicted in BN nanotubes 
[17] and BN nanoribbons p!8H2Q] . For BN nanotubes and 
nanoribbons, the transverse electric field shifts the lowest 
cr*-band downward, which is often called the nearly free 
electron state, and the nanotubes and nanoribbons can 
be metallic above a critical field [17H2Q]. The giant Stark 
effect of BN nanotubes was experimentally observed by 
Ishigami et al. [21^ . 

On the other hand, Ci et al. recently reported the 
synthesis of hybridized BN and Gr sheets using thermal 
catalytic chemical vapor deposition [22]. Since B and N 
atoms act as acceptors and donors in Gr, respectively, 
a nanoribbon made of the hybridized BCN sheet should 
have rich functionality in controlling electronic and mag- 
netic properties. Nakamura et a/.. He et al.^ and Basheer 
et al. investigated zigzag GNRs with single BN edges 
and showed that the ribbons are metallic ferrimagnetic or 



half- metallic, and that such properties are realized by the 
coexistence of extended states and edge-localized states 
[23ti25] . Kan et al. reported the half- metallic properties 
of zigzag Gr and BN hybrid nanoribbons [26 . 

In this letter, we theoretically study the electronic 
properties of the zigzag BCN nanoribbons shown in Figs. 
[l](a)-(d) using the first-principles calculations. The out- 
ermost C atoms at each edge are replaced by B or N 
atoms. Hereafter, we call this structure a B-C-N nanorib- 
bon. We show that B-C-N nanoribbons are metallic, 
since several bands cross the Fermi level. For B-C-N 
nanoribbons in a rich H2 environment, the nearly free 
electron state appears just above the Fermi level due to 
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FIG. 1. (color online) Schematic illustrations of B-C-N 
nanoribbons with N = 10 for BH,NH (a), BN,N (b), B,NH 
(c), and B,N (d). The dotted rectangles indicate unit cells. 
In Panel (e), the difference in the Gibbs free energy, AC, is 
measured from the Gibbs free energy of the B,N structure as 
a function of the chemical potential of H2 gas. 
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FIG. 2. (color online) The band structures of B-C-N nanorib- 
bons with iV = 10 for BH,NH (a), BN,N (b), B,NH (c), and 
B,N (d). The projections to p^-orbitals of B and N atoms are 
respectively indicated by the red and blue circles whose radii 
are proportional to the magnitude of the projection. 

the intrinsic giant Stark effect arising from the transverse 
electric dipole moment. We also show that the intrinsic 
giant Stark effect can be controlled by applying an elec- 
tric field. We also discuss the giant Stark effect in bilayer 
BCN nanoribbons. 

We used the projector augmented- wave method [27] 
and the local density approximation (LDA) [28 imple- 
mented in the Vienna ab-initio simulation package (vasp) 
code m [30j. As shown in Figs. [2] (a)-(d), the B-C-N 
nanoribbon is located in the x — y plane and runs along 
the X-axis. Let N be the number of zigzag lines. The cell 
size in the x-direction was chosen to be a = 2. 446 A which 
is the lattice constant of free-standing Gr optimized in 
LDA. We imposed a vacuum region about 12 A thick in 
the y and z directions. The geometry of the nanoribbons 
was fully optimized until the residual force fell below 10~^ 
eV/A. The dipole correction in the ^/-direction is used to 
exclude spurious dipole interactions between periodic im- 
ages [31 . The cutoff energy of the plane- wave basis was 
chosen to be 400 eV, and the /c-points sampling was cho- 
sen to be 16 X 1 X 1 with a Monkhorst-Pack mesh. 

First, we discuss the stability of H-terminated B-C- 
N nanoribbons. After the geometry is optimized, we 



TABLE 1. Calculated results of formation energies at zero 
temperature ^form and dipole moment in the y direction dy. 
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FIG. 3. (color online) (a) Optimized structure and electronic 
density of B-C-N nanoribbon with N = 10 for the BH,NH 
structure, (b) The same plot for the BH,N structure. The 
magnitude of projection to the pcc-orbitals of state A, which 
is indicated in Fig. [2] (a), is shown at the top of Panel (a), (c) 
The band structure of the B-C-N nanoribbon with N = 10 
for the BH,NH structure under Ey = 0.6 V/A. (d) The same 
plot for Ey — —0.6 V/A. The overlap between the cr-bands 
can be controlled by Ey^ while the 7r-bands are not sensitive 

to Ey. 

obtained perfectly flat nanoribbons independent of H 
termination. The formation energy, £^form, of a B-C-N 
nanoribbon at zero temperature is defined as 

^form = ^B-C-N " ^BN ^^^^ 

where £^b-c-n, £^bn, £^Gr, and Eyi^ are the total ener- 
gies of the B-C-N nanoribbon, BN sheet, free-standing 
Gr, and H2 molecules, respectively. Here Nq and A^h 
are the numbers of C and H atoms in the unit cell, re- 
spectively. Table 1 summarizes the formation energies 
of B-C-N nanoribbons with N = 10 for BH,NH, BN,N, 
B,NH, and B,N. We found that the BH,NH nanoribbon 
is stable at zero temperature. In this comparison, how- 
ever, the effects of environment, i.e., the temperature and 
pressure of H2 gas, are absent [32]. Stability at a finite 
temperature can be discussed by means of the Gibbs free 
energy, G, which is defined as G = £^form — A^hMHs/^, 
where /iHs is the chemical potential of H2 [32l|33]. Figure 
[1] (e) shows the difference in the Gibbs free energy, AG, 
measured from the Gibbs free energy of the B,N struc- 
ture as a function of the chemical potential of H2 gas. 
We found that the BH,NH nanoribbon is the most stable 
in a H2-rich environment and that the B,N nanoribbon 
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(a) B-C-N/B-C-N (b) B-C-N/N-C-B (c) N-C-B/B-C-N 




FIG. 4. (color online) (a)-(c) Op- 
timized structures of bilayer B-C-N 
nanoribbons with different stacking 
orders whose edges are terminated hy- 
drogen atoms. The AB stacking was 
broken, (d)-(f) The band structures 
of bilayer B-C-N nanoribbons corre- 
sponding to the structures shown in 
(a)-(c). The giant Stark effect remains 
for the B-C-N/B-C-N structure but 
disappears for the B-C-N/N-C-B and 
N-C-B/B-C-N structures. 
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is the most stable in a H2-poor environment. The BH,N 
nanoribbon becomes stablest in the intermediate envi- 
ronment, and the B,NH nanoribbon cannot be realized 
under practical conditions. 

Figure [2] shows the band structures of B-C-N nanorib- 
bons with TV = 10 for BH,NH (a), BH,N (b), B,NH (c), 
and B,N (d). The projections to p2-oi"bitals of B and N 
atoms are respectively indicated by the red and blue cir- 
cles whose radii are proportional to the magnitude of the 
projection. The highest a-band and the lowest cr*-band 
are indicated by the orange and green lines, respectively. 
The purple lines correspond to the dangling bonds of B 
or N atoms, which consist mainly of p^-orbitals. 

As shown in Fig. |2](a), the BH,NH structure is metal- 
lic, since several bands cross the Fermi level. The mag- 
nitude of projection to the p^^-orbitals of the highest a- 
band at the T point, which is labeled as A in Fig. [2] (a), is 
shown in the top Panel of Fig. [s] (a). The projections to 
the other orbitals, such as s and Py, are zero. Since the 
projection to the p^^-orbital takes a large value along the 
outermost zigzag line of the right side edge consisting of 
B atoms, the corresponding state is localized along the 
zigzag line, i.e., the current carrying edge state. The bot- 
tom of the lowest (j*-band, indicated by the green curve, 
is located just above the Fermi level. Since the effective 
mass of the lowest a* band is 1.06 mo with the electron 
mass in vacuum mo, this band corresponds to the so- 
called nearly free electron state [T7H20]. The nearly free 
electron state appears just outside of the left side edge 
consisting of N atoms similar to that induced by the gi- 
ant Stark effect in BN nanotubes and BN nanoribbons 
under an external electric field. Despite the absence of 
an electric field, the giant Stark effect seems to take place 
in the BH,NH structure as shown in Fig. [2] (a). 

The band structure of the BH,N structure is shown in 
Fig. [2] (b). We found hole pockets around the T point 
and k ^ 27r/3a. Around the F point, the highest a-band 
crosses the Fermi level, which can be regarded as the 
current carrying edge state discussed above. The purple 



curve just below the Fermi level corresponds to the dan- 
gling bond state of the N atom. We did not observe the 
nearly free electron state for the BH,N structure. Figure 
[2] (c) shows the band structure of the B,NH structure, 
where we did not observe the current carrying edge state 
but did observe the nearly free electron state. The purple 
curve just above the Fermi level corresponds to the dan- 
gling bond state of the B atom. As shown in Fig. [2] (d), 
TT and 7r*-bands of the B,N structure are quite similar to 
those calculated within tight-binding models [34]. The 
purple curves correspond to the dangling bond states of 
B and N atoms. From the band structures shown in Figs. 
|2](a)-(d), one can clearly see that the hydrogenation leads 
not only to saturation of the dangling bonds but also to 
considerable changes in the whole band structure. 

As mentioned above, the appearance of the nearly free 
electron state reminds us of the giant Stark effect in BN 
nanotubes and BN nanoribbons under an external elec- 
tric field. However, the nearly free electron state of B- 
C-N nanoribbons appears in the absence of the external 
electric field, which implies the existence of the intrinsic 
giant Stark effect due to the internal electric field, i.e., 
the dipole moment. Actually, the hydrogenated B-C-N 
nanoribbons have large dipole moments in the transverse 
direction due to the large difference in electro-negativities 
among B, C, N, and H atoms. The electro-negativities 
of B, C, N, and H atoms are 2.04, 2.55, 3.04, and 2.20, 
respectively. 

Calculated dipole moments in the ^/-direction per unit 
cell, dy^ are summarized in Table 1. The value and sign 
of dy depend strongly on the hydrogenation pattern. The 
electronic density distribution of the BH,NH and BH,N 
structures are shown by the color plot in Figs. [3] (a) and 
(b), respectively. From these figures, we found that the 
hydrogenation of the N atom leads to the formation of 
a large dipole moment in the +?/-direction along the N 
side edge of the BH,NH structure. The electric field due 
to the dipole moment is in the — ^/-direction. 

Figures [3] (c) and (d) show the effects of a transverse 
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electric field, Ey, on the band structure with = 10 
for the BH,NH structure. Applying the external electric 
field in the +^-direction to reduce the giant Stark effect, 
the nearly free electron state is shifted upward, as shown 
in Fig. [3] (c). In contrast, the application of the exter- 
nal electric field in the — ^/-direction shifts the nearly free 
electron states downward, as shown in Fig. [3](d). There- 
fore, we conclude that the nearly free electron state in 
the B-C-N nanoribbon is caused by the intrinsic giant 
Stark effect due to the internal electric field of a dipole 
moment. 

Let us move on to the bilayer B-C-N nanoribbons 
where the intrinsic giant Stark effect will depend on 
the stacking order. Figures [I] (a)-(c) show the top and 
side views of optimized structures of the bilayer B-C- 
N nanoribbons. We considered three different stacking 
orders: B-C-N/B-C-N (a), B-C-N/N-C-B (b), and N-C- 
B/B-C-N (c). The difference between the B-C-N/N-C-B 
and N-C-B/B-C-N structures is whether the outermost 
atom is N (b) or B (c). While we started the geometry 
optimization from the so-called AB-stacking structure, 
the AB-stacking was broken, as shown in the top view 
of Figs. [4] (a)-(c). The dipole moment of the B-C-N/B- 
C-N structure is dy = 0.66 eA which induces the intrin- 
sic giant Stark effect as shown in Fig. [4] (d). For the 
B-C-N/N-C-B and N-C-B/B-C-N structures, the dipole 
moments of the top and bottom layers cancel each other 
out, and the nearly free electron state around the Fermi 
level disappears as shown in Figs. [4]^ e) and (f). 

In summary, we theoretically studied the electronic 
properties of BCN nanoribbons with zigzag edges where 
the outermost C atoms are uniformly replaced by B and 
N atoms using the first-principles calculations. BCN 
nanoribbons are metallic, since several bands cross the 
Fermi level. For BCN nanoribbons in a H2-rich environ- 
ment, the nearly free electron state appears just above 
the Fermi level because of the intrinsic giant Stark effect 
due to the internal electric field of the transverse dipole 
moment. The hydrogenation of N atoms plays a decisive 
role in the appearance of the giant Stark effect in BCN 
nanoribbons. The strength of the intrinsic giant Stark 
effect can be controlled by applying an external electric 
field. With a decrease in the chemical potential of H2 gas, 
the giant Stark effect disappears due to the desorption of 
H atoms from N atoms. We also investigated the elec- 
tronic properties of bilayer BCN nanoribbons and showed 
that the appearance of the giant Stark effect depends on 
the stacking order. 
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